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We present the first in a series of microscopic studies of electrical transport through individual 
molecules with metallic contacts. We view the molecules as "heterostructures" composed of chemi- 
cally well-defined atomic groups, and analyze the device characteristics in terms of the charge and 
potential response of these atomic-groups to the perturbation induced by the metal-molecule cou- 
pling and the applied electrical field, which are modeled using a first-principles based self-consistent 
matrix Green's function (SCMGF) method. As the first example, we examine the devices formed by 
attaching two benzene-based molecular radicals-phenyl dithiol (PDT) and biphenyl dithiol (BPD)- 
symmetrically onto two semi-infinite gold electrodes through the end sulfur atoms. We find that 
both molecules acquire a fractional number of electrons with similar magnitude and spatial dis- 
tribution upon contact with the electrodes. The charge transfer creates a potential barrier at the 
metal-molecule interface that modifies significantly the frontier molecular states depending on the 
corresponding electron density distribution. For both molecules, the metal Fermi-level is found to lie 
closer to the highest-occupied-molecular-orbital (HOMO) than to the lowest-unoccupied-molecular- 
orbital (LUMO). Transmission in the HOMO-LUMO gap for both molecules is due to the metal- 
induced gap states arising from the hybridization of the metal surface states with the occupied 
molecular states. Applying a finite bias voltage leads to only minor net charge injection due to the 
symmetric device structure assumed in this work. But as current flows, the electrons within the 
molecular junction redistribute substantially, with resistivity dipoles developing in the vicinity of 
potential barriers. Only the delocalized 7r-electrons in the benzene ring can effectively screen the 
applied electric field. For the PDT molecule, the majority of the bias voltage drops at the metal- 
molecule interface. But for the BPD molecule, a significant amount of the voltage also drops in 
the molecule core. The field-induced modification of the molecular states (the static Stark effect) 
becomes significant as the bias voltage increases beyond the linear-transport region. A bias-induced 
reduction of the HOMO-LUMO gap is observed for both molecules at large bias. The Stark effect 
is found to be stronger for the BPD molecule than the PDT molecule despite the longer length 
of the former. For both molecules, the peaks in the conductance are due to electron transmission 
through the occupied rather than the unoccupied molecular states. The calculation is done at room 
temperature, and we find the temperature dependence of the current-voltage characteristics of both 
molecules is negligible. 
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I. INTRODUCTION 

Exploring the use of individual molecules as active components in electronic devices 1 has been at the forefront 
of nanoelectronics research in recent years due to the potential advantage of ultrahigh density/speed and low-cost 
device fabrication through self-assembly/self-organization processes that such molecular-scale devices may bring 2,3 . 
Numerous useful device characteristics including molecular rectifying diodes, negative differential resistance and field- 
effect transistors have been demonstrated using molecular-scale structures including small conjugated molecules 4,7 , 
single- and multi-wall carbon nanotubes 8 ' 9 and macromolecules like DNA 10 . A quantitative understanding of the 
physical mechanisms underlying the operation of such diverse molecular-scale devices has been and remains a major 
challenge in nanoelectronics research. 

The electrical characteristics of molecular devices are usually measured by sandwiching the molecules between two 
metallic electrodes 3-5,11-14 . The measured transport characteristics reflects both the nature of the metal-molecule 
interface and the properties of the molecular layer 15-25 . There are basically three electronic processes of interest 
in such metal-molecule-metal junctions: charge transfer between the metals and the molecules, the change of the 
electrostatic potential and the modification of the molecular geometry and electronic states. The nature of these 
processes under both zero and nonzero bias determines the electrical characteristics of the molecular junction. Here 
it is natural to separate the problem into device at equilibrium and device out of equilibrium since: (1) the linear 
dc-transport property of the molecular device probes the equilibrium charge distribution of the molecular junction 
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established by adsorption onto the electrodes; (2) the nonlinear transport probes the charge response of the molecular 
junction to the applied electrical field established through a finite bias voltage. A microscopic theory of molecular 
electronics will therefore need to: (1) Determine the appropriate geometry of the molecular junction; (2) determine 
the self-consistent charge transfer and the resulting lineup of the molecular levels relative to the metal Fermi-level 
and the modification of the molecular states upon formation of the metal-molecule-metal junction; (3) determine the 
molecular screening of the applied electric field, the field-induced modification of molecular states (the static Stark 
effect) and the non-equilibrium electron distribution when current is flowing; (4) determine the current/conductance- 
voltage characteristics and their correlation with the molecular and device structures. None of these issues has been 
satisfactorily solved such that theory can be used in conjunction with experiment for an unambiguous identification 
of the device operation mechanism. The purpose of our work is to elucidate the above device electronic processes and 
their dependence on the given molecular and device structures under both equilibrium and nonequilibrium conditions 
within a well-defined theoretical model. Specifically we will investigate the different aspects of molecular transport 
due to the interplay between the electronic processes at the metal-molecule interface, the electronic processses in the 
molecule core and the molecular response to the applied electrical field when current is flowing. Correspondingly, we 
focus our attention on the electron dynamics and calculate the device characteristics within the coherent transport 
regime. Our emphasis is on the conceptual understanding and the chemical trends obtained from detailed microscopic 
calculation of simple but representaive molecular and device structures. Our work can therefore be considered as a 
"minimal" quantitative microscopic model of single-molecule electronics. 

The ionic dynamics of the metal-molecule-metal junction can also be affected by the above electronic processes both 
at equilibrium and out of equilibrium, i.e., the adsorption and bias/current-induced conformation change. Neither 
aspect will be treated here. By confining ourselves to the coherent transport regime, we also neglect the effect of 
electron-vibrational coupling on transport 26 . Solving the adsorption- induced conformation change requires an accurate 
knowledge of the surface lattice structure at the atomic-scale for devices under applied voltage. This is almost never 
known in typical molecular transport measurement. Since one goal of such calculation is to provide the input nuclear 
configuration for transport calculation, our purpose can be served equally well by examining the effect that different 
adsorption and molecular geometries may have on the device characteristics, which will be treated in subsequent 
papers of the series. A rigorous solution of the bias/current-induced conformational change and the general solution 
of inelastic contribution to non-linear current due to electron-vibrational coupling is a challenging topics which requires 
solving the nonequilibrium dynamics of the coupled electron-nuclei systems 27 ' 28 since the molecular potential energy 
surface is affected by the nonequilibrium electron dynamics and electron energy is dissipated in moving the atoms. 
By focusing on the electron dynamics in the molecular junction with fixed geometry, the present work will provide a 
reference for evaluating the importance of such additional complications and provide the necessary input for further 
investigation in situations where they are indeed critical. 

In typical molecular junction measurement, conjugated molecules are attached to the metallic electrodes through 
appropriate end groups. In this paper we will consider current transport through two bencne-based molecular radicals- 
phcnylene dithiol (PDT) and biphenylenc dithiol (BPD)-adsorbcd symmetrically onto two semi-inifinite gold <111> 
electrodes through the end sulfur atoms. These structures chosen are among the simplest possible but are still 
representative of current experimental work. For the device at equilibrium, we obtain the self-consistent charge 
transfer, the adsorption-induced change in the electrostatic potential, the lineup of the molecular level relative to the 
metal Fermi-level and the modification of the individual molecular states due to the metal-molecule coupling. For the 
device out of equilibrium, we obtain the charge injection/redistribution within the metal-molecule-metal junction, the 
molecular screening of the applied field and the resulting voltage drop, the modification of the molecular states by the 
applied field, the non-equilibrium occupation of molecular orbitals, the current/conductance- voltage characteristics 
for the given contact geometry. 

The calculation we present is performed using a recently developed self-consistent matrix Green's function (SCMGF) 
method 15,16 which is based on the non-equilibrium generalization of the quasi-particle Green's function theory 29,32,33 
and uses a finite local-orbital basis set. By replacing the quasi-particle exchange-correlation self-energy 34 with 
the exchange-correlation potential within the density functional theory 35,36 , the well-established technique of self- 
consistent field theory of molecular electronic structure can then be utilized for transport calculations 16 . The real- 
space formulation of our approach allows us to provide an intuitive and coherent physical picture of the molecular 
tranport by analyzing the device elctronic processes both at the atomic-scale and on the basis of individual molecular 
orbitals. We view the molecules as comprised of chemically well-defined atomic groups and interpret their electrical 
characteristics in terms of the response of these atomic groups to the perturbation induced by the metal-molecule cou- 
pling and the applied electric field. We emphasize the insight obtained with such atomic-scale analysis and show the 
important effect of atomic-scale charge and potential inhomogeneity on device characteristics at the molecular scale, 
which may otherwise be obscured by the molecular-level analysis treating the molecules as a whole. In particular, two 
important conclusions come out of this work: (1) the adsorption-induced modification of molecular states is larger 
than the field-induced effect unless we go to large bias, so accurate modeling of the electronic processes at equilibrium 
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is critical for determining the low-bias transport characteristics; (2) the effect of both the metal-molecule coupling 
and the applied electric field in turning the individual molecular orbitals into effective conduction channels depends 
on the detailed charge and potential distribution across the metal-molecule-metal junction and may be different for 
different molecular orbitals. "Engineering" the charge and potential inhomogeneity as commonly done in quantum 
semiconductor heterostructures will be equally important at the molecular scale. Both aspects have been largely 
neglected in the past. We therefore hope the results obtained here will provide a useful guide regarding the nature 
of electron transport at the ultimate limit of device scaling and the prospect of device engineering through molecular 
design. 



II. THEORETICAL MODEL 



A. The self-consistent matrix Green's function theory 

The details of the self-consistent matrix Green's function theory have been described elsewhere 15 ' 16 , so we only give 
a brief summary here to clarify the assumptions and approximations and to show how the physical observables are 
computed. Similar but different methods based on the NEGF approach have been developed independently by other 
research groups 17,20,21 . We feel that this approach provides a natural link between quantum transport, first-principles 
electronic structure theory and qualitative molecular orbital theory 37,38 . 

Due to metallic screening in the electrodes, the charge and potential perturbations induced by molecular adsorption 
extend only over a finite region into the metal surface 39 ' 40 . We define an "extended molecule" which includes both 
the molecule and the surface atoms perturbed by molecular adsorption 41, . The surface atoms also form an adiabatic 
rcflcctionless contact with the rest of the electrodes, which can then modeled as infinite electron reservoirs commonly 
assumed for mesoscopic transport problems 16,64 . The size of the "extended molecule" is chosen such that charge 
neutrality is maintained approximately under both zero and finite biases. The central quantities in the NEGF theory 
are the retarded and correlation Green's function G r and G < 33,32 . We expand both the wavefunctions ip^ and the 
Green's functions using a finite set of local orbital functions fa , 

G^{r, ?;E)^Y. G 7 W). (1) 

i,3 

G<^(f^;E)^^2G^(E)fa(r)4,*^) (2) 

i,j 

where a is the spin index. We obtain the retarded and correlation matrix Green's function by solving the Keldysh- 
Kadanoff-Baym (KKB) equation in the matrix form 16 : 

G MM = {E + S M M - UllM ~ V Mm" ~ ^L a ( E ) ~ K'" (E)}^ 1 , (3) 

G<-'°(E) = i[G r ^(E)T L . a (E)G a ^(E)]f(E - fi L ) + z[G r ^(E)T R ., a (E)G a ^(E)]f(E - m ) (4) 
where the effect of the contact is incorporated as the self-energy operators 51 r L r R y a , 

^L^iE) — (E + Sml — Hml;o)G°ll {E + Slm — H^ M ), 



^ r R( E ) - {E + S M r - Hl IR )G RR {E + S RM - H RM ), 

L(R) — l \^L(R) ^ >L(R) 



^L(R) - ^TAm ^ r,Cr )\jR))- ( 5 ) 



Here G^, RR . — {E + Slurr) — H1 L( ^ RR ^) 1 are the surface Green's functions of the left (L) and right (R) contacts (the 
electrodes with the perturbed surface atoms removed) which can be obtained from that of the semi-infinite surface. 
H-mm represents the part of the Fock matrix contributed by the charge distributions (both nuclei and electron) in the 
"extended molecule" only, while V ext;(T (r) represents the long-range coulomb potential due to the equilibrium charge 
(ionic and electronic) distribution in the contact region, which includes the linear voltage drop due to the applied 
bias. The S are overlap matrices. The applied potential V ext changes the charge and potential in the "extended 
molecule" , requiring a self-consistent solution. 

Given the electron density or the density matrix in the "extended molecule", the calculation of the Fock matrix 
Hmm-ct is the same as that of standard molecular electronic structure calculation, which greatly facilitates the im- 
plementation of our procedure using standard molecular electronic structure codes 16 . The self-consistent calculation 



3 



proceeds by computing the input density matrix to the next iteration from the correlation matrix Green's function 
computed in the current iteration: 

The density matrix is simply the energy integration of the matrix correlation function. The integration over energy 
can be performed conveniently in the complex energy plane 16 ' 43 . Once the self-consistent calculation converges, we 
can calculate all the physical observables from the matrix retarded Green's function and the density matrix. 



B. Analyzing molecular transport 



We calculate the charge density using 



from which we can obtain the electrostatic potential in the molecular junction from the Poisson equation. We can 
also calculate charges associated with each atom using Becke's atomic-partition scheme 44 , 



N 



J drp(r) = £ / drWi{f)p{f) = ^ N t (8) 



Here the atomic weight function Wi{r), which satisfies J2 i Wi(r) — 1 everywhere in space, is determined such that it 
is centered on the atom i and is non- negligible only in a region close to its atomic center 44 . 
The local density of states (LDOS) gives the energy-resolved charge density distribution: 

n CT (r, E) = -- lim £ Imag[G%°(E + iS)]M^{f), (9) 

The spatial integration of LDOS gives the density of states, 

n a (E)= ( drn a {r,E) = -- lim Tr{Imag[G R ' a {E + iS)]S}, (10) 

J 7T <5-^0+ 

To identify the contribution of the individual molecular orbitals to the total density of states, we project it onto the 
basis of molecular orbtias. This is done by transforming G R into the basis of the molecular orbitals obtained by 
diagonalizing the (self-consistent) Fock matrix corresponding to the molecule. The imaginary part of the diagonal 
element of the transformed G R gives the projected density of states (PDOS) of the corresponding molecular orbitals. 
The peak position of the PDOS characterizes the perturbed molecular energy level, while the broadness reflects the 
coupling strength of the molecular oribitals with the metallic surface states. 

A critical question in molecular transport is how the electron occupation of the molecular orbitals changes as the 
molecule is driven out of equilibrium by a finite bias voltages. The information is contained in the non-equilibrium 
density matrix p°- (Eq. 6). Transforming to the basis of molecular orbitals, the diagonal elements of the density 
matrix give the the non-equilibrium electron occupation of the corresponding molecular levels. 

A general formula for the current through a mcsocopic system with arbitrary interaction in contact with two 
non-interacting electrodes is: 

I= \\ dEj2Tr{[(Tl-T° R )(E,V) iG<^(E,V)\ 
+ [f(E - (i L )ri(E, V) - f(E - ii R )r> R (E, V)]A°(E, V)}} (11) 

where A a = i(G r;<T — G a;<T ) is the spectral function. Since the only scattering mechanism in the coherent transport 
regime is that by the contacts (introducing additional scattering mechanisms within the molecule will lead to additional 
terms in self-energy and current), we arrive at the familiar Landauer-type current formula: 

I= lj dEj2T°(E,V)[f(E-» L )-f(E- m )} (12) 
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where the transmission probability T a through the molecule 45 ' 46 ' 16 is obtained from: 

T a (E, V) = Tr[i a {E, V)] = Tr[Y a L {E, V)G R '' a {E, V)T a R {E, V)[G R ' a ]\E, V)], (13) 

Note that equations (12-13) hold for both orthogonal and non- orthogonal basis functions. To identify the contribution 
of individual molecular orbitals to the transmission, we transform the transmission matrix t a (E, V) into the the 
basis of the molecular orbitals whose diagonal elements give the transmission probability through the corresponding 
molecular orbitals. 

Combined with the spatially resolved LDOS and charge density, the total DOS, the transmission coefficient and 
their projection onto individual molecular orbitals provide a set of useful qualitative analysis tools to establish the 
connection between the molecular electronic structure and the transport characteristics of the metal-molecule-metal 
junction, similar to the use of qualitative molecular orbital theory in quantum chemistry 38 . A caveat of projecting 
the physical observable onto the individual molecular orbitals is that phase interference information is lost. Only 
the summation over all molecular orbitals (conserved by the matrix transformation), not the individual components, 
corresponds to the quantum mechanical physical observable. The situation is similar to what we met when we try 
to decompose the electron density of the molecule into the contribution of individual atoms through a Mullikcn- 
type population analysis or other charge partition schemes as we used here. Only the total density but not the 
atom-partitioned density corresponds to a physical observable. 

Within the coherent transport model, the only temperature effect is from the two electrode Fermi-Dirac distri- 
butions. We can separate the current into two components, the "tunneling" component I tun and the "thermionic 
emission" component Ith as follows, 

p run ft*L r+oo 

I = Itun+Ith = -r[ +(/ +/ )}dET(E,V)[f(E - »l) - f(E - » R )} (14) 

J HL J -OO J fJ, R 



C. Device model 

We choose the electrostatic potential in the middle of the (empty) bimetallic junction as the energy reference, then 
the equilibrium Fermi-level Ef is fixed at the negative of the metal work function —5.31(61^) for single-crystal gold 
electrodes. The electrochemical potential of the two electrodes is fixed by the applied bias voltage V at hl — Ef — eV/2 
and fiR = Ef + eV/2. Note the Fermi-level positions are fixed by the bimetallic junction alone without the molecular 
insertion. The voltage drop across the molecular junction is determined by the molecular response to the metal- 
molecule coupling and the applied bias voltages. The bias polarity is chosen such that for positive bias the electron 
is injected from the right electrode. 

In this work, the electronic structure of the "extended molecule" is described using the Becke-Perdew- 
Wang(BPW91) 48,49 parameterization of the spin-polarized density-functional theory (SDFT) 35,50 within the 
Generalized-Gradient Approximation (GGA) 49 . We also replace the atomic core by an appropriate pseudopoten- 
tial 51 with the corresponding optimized Gaussian basis sets 52,53 . The geometry of the adsorbed molecule is taken to 
be the same as the singlet geometry of the free molecule optimized at the BPW91/6 — 31G* level (we asuume the bare 
molecule to be the molecular biradical with the end H atoms removed) . The adsorption geometry is chosen such that 
the end atoms sit in front of the center of the triangular pad of the three gold atoms on the Au<lll> surface (the end 
sulfur-suface distance is 1.9 A). Six nearest-neighbor gold atoms on each metal surface are included into the "extended 
molecule" . Within the range of second-nearest-neighbor coupling, there are 12 metal atoms in the first surface layer 
and 14 metal atoms on the second surface layers on each side coupled to the "extended molecule" . Only the blocks 
of the surface Green's function corresponding to these atoms are needed, which in turn can be calculated from the 
surface Green's function of the semi-infinite metal using tight-binding method parametrized by fitting accurate bulk 
band structure calculations 15,54 . The results and conclusions given in the following are not affected significantly by 
small changes in the molecular and adsorption geometry. The structures of the molecule junction are shown in Figs. 
(1) and (2). The calculation is performed using a modified version of GAUSSIAN98 15,55 . 

III. DEVICE AT EQUILIBRIUM 

The problem at equilibrium is a generalization of the familiar chemisorption problem in surface science since two 
metallic surfaces are involved. For a molecular orbital to be turned into an effective conduction channel of the metal- 
molecule- metal junction (within the coherent transport regime), it needs to couple well to both electrode states. 
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Its energy must lie within the energy window determined by the bias voltage and the thermal broadening through 
the electrodes Fermi-distribution. Both aspects are affected by the molecular adsorption, which may well be more 
significant than the subsequent application of the bias voltage, as detailed in subsequent sections. 

A. The electronic structure of the molecules and the identification of chemical groups 

The transport property of the molecular junction is determined mainly by the hybridization of the surface metal 
states with the frontier molecular orbitals, i.e., the molecular states whose energies lie closest to the metal Fermi- 
level. So we start by illustrating in Fig. (3) and Fig. (4) the charge distribution of the HOMO-l,HOMO,LUMO and 
LUMO+1 states of the two isolated molecules (with the end H removed). For both molecules, the HOMO-1 states are 
localized on the end sulfur atoms, the HOMO states are delocalized over entire conjugated backbone while the LUMO 
and the LUMO+1 states are benzene based. The LUMO state of BPD also has a hnite weight on the end sulfur 
atoms. Both molecules can be viewed as "heterostructures" composed of the end sulfur atoms and the benzene rings. 
Due to the non-coplanar geometry of the two benzenes in BPD molecule (the torsion angle is 37°), the orbital overlap 
between the corresponding it electrons is weak. Our analysis of the electrical transport through the two molecules 
will be based on the electrical response of these chemical units. 

B. Charge transfer and electrostatic potential change in the molecular junction 

The adsoption-induccd charge transfer across the metal-molecule interface is the central quantity for the device at 
equilibrium since the linear-response transport probes the equilibrium charge distribution of the molecular junction. 
The perturbation introduced by the metal-molecule coupling is largely a localized interaction. The corresponding 
charge-transfer process is determined by the interfacial chemistry and can be understood qualitatively from the local 
bonding analysis across the metal-molecule interface as in the chemisorption problem 15,37 . 

Due to the identical end group and bonding configuration across the metal-molecule interface, both the magnitude 
and the spatial distribution of the transferred charge for the two molecules are quite similar. This is illustrated in Figs. 
(5-6), where we plotted the difference between the self-consistent charge distribution in the gold-moleculc-gold junction 
and the charge distribution in the isolated molecule plus the charge distribution of the isolated bimetallic contact. 
For clarity and to aid visualization, we have plotted both the regions where charge accumulation and depletion occur 
and the spatial distribution of the transferred charge as a function of position in the X-Y plane (defined by the left 
benzene ring). The charge transfer process involves mainly the end sulfur atoms and neighbor carbon atoms and 
decays rapidly as we move away from the metal-molecule interface. This is clearly seen for the BPD molecule where 
the charge perturbation induced by the two metal-molecule contacts decays into the interior of the molecule without 
effectively interfering with each other (Fig. 6). Note this is not necessarily due to the non-coplanar geometry. The 
decay of charge perturbation is already obvious for the PDT while for BPD it becomes negligible before reaching 
the inter-benzene bonding region. The number of electrons increases in the sulfur-gold bonding region due to the 
rehybridization of the sulfur P z orbital in forming gold-sulfur bond and the transfer of charge from gold atom to 
sulfur atom due to electronegativity difference. Electron density also increases in the neighbor carbon P z orbital since 
the HOMO level has large weights on both the sulfur and carbon P z orbitals. The electron density decreases in the 
sulfur-carbon c-bonding region, the sulfur P x orbitals and also the gold s orbitals. This is because the formation of 
surface gold-sulfur bond involves charges originally residing in the sulfur P x and gold s orbitals thus weakening the 
bonding between sulfur and carbon. The direction of charge transfer is from the gold electrodes to the molecules. 

Accompanying the charge transfer across the gold-molecule interface, the electrostatic potential also changes. The 
difference of the electrostatic potential in the junction and the electronic potential in the isolated molecule plus that 
in the isolated bimetallic contact gives the potential perturbation due to the formation of the contact. This is plotted 
in Figs. (7) and (8). The net transfer of electrons into the molecule increases the electrostatic potential inside the 
molecule creating a potential barrier between the metal surface and the end sulfur atom. There is also a potential well 
between the sulfur atom and the benzene ring because of the decrease of electron density in the sulfur-carbon bonding 
region. Although the pattern of charge transfer at the metal-molecule interface is similar for the two molecules, the 
long-range electrostatic potential perturbation is quite different. For the BPD molecule, the electrostatic potential 
profile inside the the two benzene rings becomes quite complicated which creates additional barrier for the electron 
motion within the molecule core. Due to the non-coplanar geometry, the barrier for injection from the right metal 
into the right half of the molecule is larger than that from the left metal to the left half (the right benzene ring 
gives a slightly smaller orbital overlap with the right electrode). Since the probability of electron tunneling through a 
potential barrier is sensitive to its shape, the charge-transfer induced potential change will have a profound influence 
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on the electron transmission through the metal-molecule-metal junction. The presence of the potential barrier at the 
interface (and also in the molecule core for BPD) will affect the charge redistribution within the molecule when an 
additional field is applied, since it impedes the flow of electrons across the metal-molecule-metal junction. This will 
be discussed in the next section. 



C. Contact-induced modification of molecular states, band lineup and conductance of the 

metal-molecule-metal junction 

The coupling between the gold electrodes and the molecule (reflected in the self-energy operator of the electrodes) 
and the induced potential change in the the molecule (reflected in the self-consistent molecular Fock matrix) modify 
both the charge distribution and the energy of the molecular states relative to the metal Fermi-level. They also broaden 
the discrete molecular electronic spectrum into a quasi-continuous one. These effects can be analyzed through the 
local density of states, the transmission coefficient and their projection onto the individual molecular orbitals. 

The transmission versus energy (T-E) and the projected DOS (PDOS) corresponding to the five frontier molecular 
states of the two molecules (LUMO+1, LUMO, HOMO, HOMO-1 and HOMO-2) closest to the Fermi-level for spin- 
up electrons are shown in Figs. (9) and (10). The molecules studied here contain an even number of electrons 
and the gold electrode is non-magnetic, so the spin-up and spin-down channels show identical characteristics. For 
other molecules in contact with ferromagnetic electrodes, the spin-resolved analysis is essential for understanding the 
transport characteristics. 

From both the transmission versus energy curve and the projected DOS curve, we find the Fermi-level lies closer to 
the HOMO state than to the LUMO state for both molecules. Comparing the peak position in the PDOS plot with 
the energy levels of the isolated molecule, we find that the contact with the metallic electrodes significantly shifts 
the energy levels of the frontier molecular states. The change is found to be larger for the occupied states than for 
the unoccupied states since the unoccupied states are benzene based (which also leads to sharper peaks in the PDOS 
plot). 

The contact with the electrodes also modifies the charge distribution of the molecular states. This is illustrated 
by examining the LDOS at energies corresponding to the peak positions of the PDOS of the HOMO and LUMO for 
both molecules in Figs. (11) and (12). We show both the shape and the spatial distribution of the surface-perturbed 
molecular states for PDT. The spatial distribution is obtained by integrating the 3-D LDOS with respect to the z-axis 
and plotted as a function of position in the xy-plane. The perturbation of the molecular states due to the coupling 
to the electrodes is different for the two molecules. In general, the surface-induced change in the charge distribution 
associated with the molecular states correlates closely with the change in the electrostatic potential. For the PDT 
molecule, the large increase of the electrostatic potential in the middle of the benzene rings pushes the electrons 
located on the carbon atoms to the peripheral hydrogen atoms for both the HOMO and the LUMO states (Fig. 11). 
The LUMO level remains highly localized giving rise to sharp peak in both PDOS and T-E plots. For the BPD 
molecule, electrons also move from the right benzene to the left benzene for both the HOMO and LUMO states due 
to the smaller potential arising from the difference in local geometries. Similar analysis applies to other molecular 
states. 

The peak positions in the PDOS plots of the HOMO and LUMO states align closely with the two peak positions 
in the T-E plot around the Fermi-level, corresponding to the onset of resonant transmission. By projecting the 
total transmission coefficient onto the individual molecular orbitals, we find the first transmission peak above Ep (at 
E = — 2.45(eU)) arises mainly from the nearly degenerate LUMO and LUMO+1 states. The first transmission peak 
below Ep (at E = — 6.5(eU)) instead arises mainly from the HOMO state. In addition to the transmission peak at 
the HOMO and LUMO states, there is also a transmission peak in the HOMO-LUMO gap of the PDT molecule at 
energy E — — 3.6(eU) (Fig. 9). This peak corresponds to transmission through metal-induced gap states (MIGS) due 
to the hybridization of metal surface states with the occupied molecular states, which is significant for such a short 
molecule as PDT. The projected transmission analysis shows contribution coming mainly from the HOMO and also 
other occupied molecular states. The nature of the MIGS state is most clearly seen in the corresponding LDOS plot 
(Fig. 13), which shows similar charge distribution to the HOMO state. This is the case for both the LDOS and the 
transmission coefficient throughout the HOMO-LUMO gap region including that at the Fermi-level. For such a short 
molecule as PDT, the transmission through the metal-induced gap states in the HOMO-LUMO gap is significant. A 
similar situatuation occurs also for the BPD molecule except that the transmission through the HOMO-LUMO gap of 
the BPD molecule is much reduced (Fig. 10). The molecule is longer and the orbital overlap with the right electrode 
states is weaker, so the transmission probability in the gap is much smaller. 

The magnitude of the transmission coefficient at the Fermi-level determines the zero-bias conductance at low 
temperature. We find a conductance of 4.8 (f/,S) and 1.4(/zS) for the PDT and BPD molecules. Since the length of 
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the molecules are 6.4(A) and 10.7(A) respectively, the resistance is not proportional to the conjugation length, as 
expected for tunneling transport. 

IV. DEVICE OUT OF EQUILIBRIUM 
A. Current-voltage and conductance-voltage characteristics 

The current-voltage (I-V) and the differential conductance- voltage (dl/dV-V or G-V) characteristics are calculated 
for the two molecules using the method described in Sec. II in the bias range from — 4(V) to 4(V) and plotted in Figs. 14 
and 15. Due to the symmetric device structure, the current-voltage and conductance- voltage characteristics for both 
molecules are nearly symmetric with respect to bias polarity. In the low bias regime, current changes approximately 
linearly with the bias voltage for both molecules, corresponding to tunneling transport in the HOMO-LUMO gap. Both 
tunneling and thermionic emission contribution to the total current are shown. As expected, for tunneling transport 
through a large barrier, the thermionic emission contribution to the current is negligible even at low bias (shown in 
the insets). Since the variation of the thermionic emission contribution with bias voltage is small (in the coherent 
transport regime), this further reduces its contribution at high bias. Since this is the only temperature dependence 
in the coherent transport model, we expect the temperature dependence of the coherent current transport through 
the two molecules to be weak (ignoring any disorder effects). Within the coherent transport model, the thermionic 
emission and correspondingly the temperature dependence of the current transport can be important only if the the 
barrier for electron transmission is small. 

For comparison, we have also plotted the conductance-voltage characteristics obtained using the equilibrium trans- 
mission coefficient, .i.e., replacing T(E, V) in Eq. 12 with T{E 1 V = 0). Since the gold surface density of states 
are approximately symmetric with respect to the Fermi-level (the gold surface band around the Fermi-level is due 
primarily to the sp electrons), the difference between the G-V characteristics thus obtained and the self-consistent 
G-V characteristics reflects the effect of the bias-induced modification of molecular states. 

For both molecules, only the low-bias conductance-voltage characteristics (before reaching the first conductance 
peak) can be reasonably well reproduced by the equilibrium transmission characteristics. The deviation in both the 
magnitude and the peak position of the conductance becomes significant at large bias, there are also more conductance 
peaks than would be obtained from the equilibrium transmission characteristics. So any attempt to predicting the 
nonlinear transport characteristics from the equilibrium transmission characteristics combined with an assumption 
about the voltage drop will lead to significant error. The effect of the bias-induced modification of molecular states 
is also obvious from looking at the shift of the frontier molecular levels by the applied voltage, as shown in Fig. (16). 
The molecular levels plotted are obtained by diagonalizing the molecular part of the self-consistent Fock matrices at 
each bias voltage. The molecular levels are nearly constant at low voltages, but begin to shift before reaching the 
first conductance peak. The voltage at which the shift occurs corresponds to the voltage where major deviation in 
the G-V characteristics occurs. 

An important question in current transport through molecules is which molecular states are responsible for the 
observed conductance peak. From the position of the plotted molecular levels, we expect that for both molecules, 
the peaks in the conductance are due to the occupied molecular levels, in disagreement with previous calculation 
where a jellium model of the electrode is used 18 ' 19 . We believe this is because the jellium model underestimates the 
metal work function which effectively pushes the HOMO level down relative to the metal Fermi-level 15 . The shifts in 
the individual molecular levels are not identical in either direction or magnitude, so a rigid shift in energy levels as 
occuring in a planar metal-semiconductor contact does not apply here. Instead the shift of the individual molecular 
states will depend on the detailed charge and potential distributions. For the two molecules considered here, there is 
an effective reduction in the HOMO-LUMO gap as bias increases (Fig. 16). 

B. Charge response and voltage drop 

The charge response of the molecule to the applied electrical field is reflected in both the net charge injection into the 
molecule and the charge redistribution inside the molecule. For PDT molecule, we find there is only slight additional 
charge injection into the molecule as we apply the bias voltage. For BPD molecules, the net charge injection becomes 
important only when we go to high bias (> 2.5(eV)). The charge injection under nonzero bias voltage is determined by 
the balance of charge injection and extraction at the source-molecule and drain-molecule contacts. Since for the PDT 
molecule the contact geometries are identical and the gold surface density of states are approximately symmetric with 
respect to the Fermi-levels, little net charge accumulation will be induced by the applied bias voltages. For the BPD 
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molecule, the two rings are non-coplanar, but due to the largely localized nature of the charge injection process, the 
difference in the two contact becomes non-negligible only at high bias voltages. In addition to the total net charge in 
the molecule, we can also calculate the noncquilibrium occupation of the individual orbitals from the nonequilibrium 
density matrix (Eq. 6) as a functon of the bias voltage (Fig. 17). We find the change in the electron occupation is 
gradual at the voltages corresponding to the conductance peak, where molecular level moves past the fermi level of 
the electrodes. The electron occupation of the molecular orbital is fractional due to the broadening of the molecular 
level by the coupling to the contact. 

Although the net charge injection due to applied voltage is negligible, there can be significant charge redistribution 
within the molecule, hiscalls for closer investigation into the spatial distribution of the charges and potential distri- 
bution within the molecule as a function of applied voltage 63 . Samples of the spatial distributions of the transferred 
charge and the electrostatic potential drop for both molecules at bias voltage of 3.0(V) are shown in Fig. 18. Similar 
charge and potential distributions are obtained for all bias voltages. The spatial distribution of charge transfer is 
obtained by integrating the difference in electron density at finite and zero biases along the z-axis and plotted as a 
function of position in the xy-plane (defined by the left benzene ring). The potential drop is obtained by evaluat- 
ing the difference between the electrostatic potential at finite and zero biases, which obeys the boundary condition 
of approaching +V/2 (— V/2) at the left (right) electrode. The molecules exhibit different conductance at different 
voltages, but the calculated charge and potential distributions show similar patterns and do not depend on whether 
the molecules are in low or high conductance state. The reason is as follows: The charge and potential response of the 
molecular junction is determined by the total electron population. But to be in high conductance state, one molecular 
orbital needs to align with one of the metal Fermi- levels. As it moves away from alignment, the change in the electron 
population is gradual (Fig. 17), in contrast with the change in conductance. Given the net charge injected into the 
molecule, the determination of the electrostatic charge and potential response of the molecule in contact with the two 
electrodes will be equivalent to that of an isolated charged molecule under the same boundary conditions of charge 
and potential. 

The nature of the charge redistribution can be understood readily by partitioning the molecules into the atomic- 
groups, i.e., the benzene rings and the end sulfur atoms. As bias voltage increases, electrons move from the source-side 
carbon and sulfur to the drain-side carbon and sulfur. This is due to fact that the tt electrons in the carbon Pz orbitals 
can move freely under the applied electric field. This flow of tt electrons is impeded at the molecule-drain contact 
because the the presence of potential barrier there, which inhibits the charge flow into the drain electrode. Similar 
considerations apply to the molecule-source contact. This leads to the creation of two resistivity dipoles at the 
metal-molecule interface, i.e., the accumulation of electrons on the injecting side and the depletion of electrons on 
the extracting side of the potential barrier. For the BPD molecule, where the two benzene benzene rings are non- 
coplanar, the charge response of the two benzene rings show identical behavior and are similar to the benzene ring in 
the PDT molecule. The weak orbital overlap between the two benzene rings disrupts the flow of tt electrons across 
the inter-benzene bonding region, which creates another potential barrier with corresponding resistivity dipole and 
partially insulates the two benzene rings from each other. 

The spatial variation of the current-induced electron redistribution could also have a significant effect on the 
structural stability of the molecule under high bias. As shown in Fig. 18, applying bias voltage increases the electron 
density in the sulfur-carbon bond at the source-molecule contact but decreases the electron density in the sulfur- 
carbon bond at the drain-molecule contact. According to Hellman-Feynman theorem 66 , this would lead to stronger 
attractive forces and consequently shorter sulfur-carbon bond at the source-molecule contact but weaker forces and 
longer sulfur-carbon bond at the drain-molecule contact. For the BPD molecule, there is also a shift of electron 
density in the carbon-carbon bond connecting the two benzene rings which may affect the inter-ring spacing. Since 
the magnitude of the electron redistribution increases with the bias voltage, we expect this bias-induced modification 
of structural change will have stronger effect at high-bias. The problem of current-induced conformational change is 
the molecular analogue of the more familiar electromigration problem in metallic systems . Although this problem 
has been treated recently by several groups 19 ' 68 , we feel that much needs to be done before a satisfactory theory 
emerges. 

The resistivity dipole is a well known concept in mesoscopic electron tranport 64 and is common in inhomogenous 
transport media where their presence in the vicinity of local scattering center helps to overcome the barrier for trans- 
port and ensure current continuity. This can lead to a nonlinear transport effect due to the strong spatial variations 
in local carrier density and transport field. Such nonlinear transport characteristics induced by the device charge and 
potential inhomogeneity are a hallmark of band "engineering" through mesoscopic semiconductor heterostructures 65 . 
Molecular junctions are an analogue of the mesoscopic "heterostructures", since the choice of the component func- 
tional and structural groups allows the possibility of "engineering" charge and potential inhomogeneity within a single 
molecule, suggesting the possibility of device "engineering" through molecular design and permitting the extension of 
device concepts to the molecular scale. This effect can be further elucidated by examining the potential response of 
the two molecules. 
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The mobile 7r-electrons of the benzene ring effectively screen the applied field, leading to rather flat electrostatic 
potential drop across the benzene rings. The electrostatic potential drops rapidly in the molecule-electrode contact 
region, and becomes flat again as it approaches the boundary of the electrode where strong screening of the applied 
electric field occurs. For the PDT molecule, the majority of the voltage therefore drops across the molecule-electrode 
contact region 17,63 . Since the unoccupied molecular states of the PDT molecule are located mainly on the benzene 
ring where the screened electric field is small, their energy levels don't change much as bias voltage increases. The 
HOMO-1 state shows stronger voltage dependence than the HOMO state because it is end-sulfur based where the 
potential variation is the strongest. For the BPD molecule, although the ir electrons of each benzene ring can screen 
the applied field effectively, the potential barrier between the two benzene rings prevents the electron from flowing 
freely across, leading to significant voltage drop across the carbon framework. As a result, there exists strong spatial 
variation of the electrostatic potential across the molecule core in addition to the metal-molecule contact. Since the 
four frontier molecular states of the BPD molecule are either sulfur-based or have charge distributed across the entire 
molecule core, they all show clear voltage dependence. An empirical model of rigid shift of the molecular level relative 
to the Fermi-levels of the two electrodes has often been assumed 11 , with the proportion of the voltage drop being 
associated with the contact geometry. It is clear that this model can only be utilized as a crude check of the device 
characteristics, since it neglects the possible voltage drop across the molecules and the different effects this may have 
on different molecular states. 



C. Bias-induced modification of molecular states and transmission characteristics 

The molecular charge and potential response to the applied electrical field will affect the transport characteristics in 
two ways: (1) it shifts the molecular level relative to the Fermi-levels of the two electrodes; (2) it modifies the charge 
distribution of the molecular states and therefore their capability for carrying current. The bias- induced modification 
of molecular states can be analyzed by examining the local density of states and its projection onto the individual 
molecular orbitals at finite bias voltages. 

A snapshot of the bias dependence of the transmission characteristics and the projected DOS corresponding to the 
five frontier molecular orbitals is shown for the PDT molecule at bias voltages of 1.6(F), 3. 0(V), 4.0(F) (Fig. 19) and 
for the BPD molecule at bias voltages of 1.4(F), 3. 0(F), 4.0(F) (Fig. 21) to illustrate their voltage dependence. We 
have also shown for PDT molecule the LDOS at energies corresponding to the peak position in the PDOS plots of 
the HOMO and LUMO states at bias voltage of 3.8(F) (Fig. 20). 

For the PDT molecule changing from zero bias (Fig. 9) to V = 1.6(F) (Fig. 19), both the peak positions in the 
PDOS and the electron distribution associated with the LUMO states don't change much. The HOMO energy level 
also doesn't change much, but there is a shift in charge distribution from the right end sulfur atom to the right carbon 
atom (not shown here). The energy of the HOMO-1 state increases with bias, leading to a decrease in the energy 
spacing between the HOMO and HOMO-1 states. At 1.6(F), the transmission versus energy of the PDT molecules 
reaches a peak at energy corresponding to the drain Fcrmi-lcvcl at Ep — eF/2, giving rise to the first peak in the 
conductance. As we increase the bias voltage further from F = 1.6(F) to 3.0(F) (Fig. 19), the energy levels of 
the HOMO and HOMO-1 states shift gradually toward the equilibrium Fermi-level Ep, increasing the transmission 
coefficient there. As we further increase the bias to V = 3.8(F), the charge distribution (the LDOS) at energies 
corresponding to the peak position of the LUMO state also changes, developing large weight on the end sulfur Pz and 
the peripheral hydrogen atoms (Fig. 20). Because the HOMO and HOMO-1 states move up to the metal Fermi-level 
as the bias increases, the current increases more gradually than would be obtained if the bias-induced modification 
of molecular states is neglected (Fig. 14). 

For the BPD molecule, the first peak in conductance is reached at V = 1.4(F) (Fig. 15) and Fig. 21). From 
zero bias to 1.4(F), the energies of the LUMO and HOMO-1 states decrease and there is a shift of charge from the 
right benzene to the left benzene ring. The energy of the HOMO states increases, and the shift of charge distribution 
is from the left benzene ring to the right benzene ring (not shown here). Compared to the PDT molecule, the 
bias-induced modification is stronger despite the longer molecule length, due to the stronger spatial variation of the 
potential profile. As we increase the bias to V = 3.0(F) (Fig. 21), the energy of the LUMO continues to decrease, 
but the energies of both the HOMO and HOMO-1 states increase toward the metal Fermi-level Ep. The transmission 
coefficient at both energies decreases significantly since the electrons are now more localized on the right ring. A 
notable change is that a second peak in the PDOS of the HOMO state develops at E = — 7.0(eF). Examining the 
corresponding LDOS shows similar charge distribution to that of the HOMO-2 state with energy of — 7.15(eF) which 
has charge distribution on both benzene rings (not shown here) and correspondingly large transmission probability. 
Increasing the bias to 4.0(F) (Fig. 21) further decreases the energy of the LUMO state and increase the energies of the 
HOMO and HOMO-1 states. At F = 3.8(F), the electrons shift from the right benzne ring back to the left for both 
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the HOMO and HOMO-1 states (not shown) increasing the transmission coefficients there. Since the electrochemical 
potential of the left (right) electrode approaches alignment with the HOMO-2 and also the LUMO state, there is a 
rapid increase in the conductance as we increases the bias further toward V = 4.0(V) (Fig. 15). 

V. DISCUSSION AND CONCLUSION 
A. Electron transport or hole transport? 

A common practice in the literature on molecular electronics is to characterize the molecular transport as "electron 
transport" if the conduction is mediated by tunneling through the unoccupied molecuar states, or as "hole transport" 
if the conduction is mediated by tunneling through the occupied molecular states, following the terminology commonly 
used in semiconductor/organic device and electron transfer research. Here it is important to recognize their differences. 

For bulk semiconductor devices, the concepts of "electron" or "hole" are associated with introducing shallow- 
impurity (dopant) atoms into the ideal lattice structure, which introduces electron intos the (unoccupied) conductand 
band or removes electrons from the (filled) valence band. The "electron" or "hole" thus introduced are mobile charge 
carriers with delocalized wavefunctions. The system as a whole remains charge neutral, since the charges associated 
with these carriers are compensated by the impurity ions left behind. Therefore, the type of charge transport is an 
equilibrium property depending only on the type and amount of dopant atoms introduced. By contrast, for organic 
devices, the system is often not doped. Charge transport occurs by injecting charge carriers into the system from 
the electrodes. The type of transport depends on the charging state of the molecule and the nonequilibrium injection 
of charge carriers into the neutral system. For both cases, the type of transport is determined by the change in the 
occupation of electron states through doping or contact injection, indpendent of the nature of charge transport itself 
(band transport or polaron hopping). 

The situation for coherent molecular transport is quite different. For the symmetric molecular devices considered 
here, both the total number of electrons (the charging state of the molecule) and the occupation of the individual 
molecular orbitals change little with the applied bias. Although a fractional amount of charge is transfered from the 
gold to the molecule upon electrode contact, the transfered charge is localized in the interfacial region and characterizes 
the changes in the interfacial bond. The change in the occupation of molecular states is fractional and gradual, cither 
upon contact to the electrodes at equilibrium or upon application of a nonzero bias voltage out of equilibrium (Fig. 
(17)). This is due to the quantum mechanical nature of the coherent tunneling transport, where we cannot characterize 
the tunneling electron as being physically injected into the molecule and subsequently extracted out. Therefore, it 
is more appropriate to characterize molecular transport through the resonant molecular states without associating it 
with specific transport types. 

B. Single molecule versus molecular monolayer 

The focus of the present series of work is on current transport through a single molecule in contact with two metallic 
electrodes. This corresponds closely to the molecular transport measurements using atomic-size break junctions 4,69-72 , 
where either one or several molecules are probed. However, many molecular transport experiments are performed 
on molecular monolayers where thousands or tens of thousands of molecules are probed by the contacts to the 
electrodes 3 ' 6,7 ' 73 . Such experiments with monolayer configuration present a quite different situation from the single 
molecule configuration considered here. Besides the additional complexity of inter-molecular interactions within the 
monolayer, the most important difference lies in the interface electrostatics. Note that the same boundary condition 
for the electrostatic potential across the molecular junction applies to both the single molecule configuration and the 
monolayer configuration: deep inside the electrodes they approach the bulk value, which can be shifted rigidly with 
respect to each other by the applied bias voltage. For the single molecule configuration, the transfer of charge across 
the metal-molecule interface is confined in a small region, whose contribution to the electrostic potential decays to 
zero in regions far away from the molecule. But for the monolayer configuration, the electrostatic potential in regions 
far away from the molecule is the superposition of contributions from the transfered charge on a large number of 
molecules. To satisfy the same boundary condition of the electrostatic potential, the charge transfer per molecule 
in a molecular monolayer can be orders of magnitude smaller than that in the single molecule considered here 74,75 . 
This is the situation often met in organic electronics, where a monolayer of self-assembled molecules is used to modify 
the work function of the metallic contacts 61 ' 62 . Similar problems have also been considered in the electron transport 
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through metal-carbon nanotube interfaces ' . Here it is important to recognize the different physics of the metal- 
molecule interface in single-molecule devices and monolayer devices, since it may have profound effects on the device 
functionalities achievable using molecular materials. 



C. Conclusion 



We have presented a first-principles based microscopic study of current transport through individual molecules. The 
real-space formulation allows us to establish a clear connection between the transport characteristics and the molecular 
electronic structure perturbed by the metal-molecule coupling and the applied electric field. By separating the problem 
into equilibrium and non-equilibrium situations, we identify the critical electronic processes for understanding the 
linear and non- linear transport characteristics. At equilibrium, the critical problem is the charge transfer process 
upon formation of the metal-molecule-metal contact, which modifies the molecular states and determines the energy 
level lineup relative to the metal Fermi-level. This is mainly an interface-related process and can be controlled by 
controlling the contact. Out of equilibrium, the central problem is the molecular charge response and the consequent 
molecular screening of the applied electric field, which depends on both the molecule core and the nature of the metal- 
molecule contact and can be understood by viewing the molecules as "heterostructures" of chemically well-defined 
local groups and analyzing their electrical response to the applied electrical field 78 
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example, the calculations reported in this work were performed using a quite restricted basis set for the metal atoms on 
the surface of the semi-infinite substrate. The metal atoms belonging to atomic-scale structures on top of the substrate 
and the molecule itself are treated using much better basis sets as in standard molecular calculations. It is clear that the 
description of the metallic substrate is inadequate. And obtaining a better balance in the basis set decription between the 
substrate atoms and the molecule is a difficult and important issue for further improving the accuracy of the modeling of 
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FIG. 1. Atomic geometry of the gold-PDT-gold junction. Six gold atoms closest to the end sulfur atoms on each electrode 
are included into the "extended molecule". 

FIG. 2. Atomic geometry of the gold-BPD-gold junction. Six gold atoms closest to the end sulfur atoms on each electrode 
are included into the "extended molecule". 



FIG. 3. Orbital shape of the HOMO-1, HOMO, LUMO and LUMO+1 states of PDT 
FIG. 4. Orbital shape of the HOMO-1, HOMO, LUMO and LUMO+1 states of BPD 

FIG. 5. Charge transfer upon the formation of the gold-PDT-gold contact. The top and the middle figure show the isosurface 
plot of the region where electron increases and decreases respectively. The bottom figure shows the spatial distribution of the 
transfered electrons as a function of position in the xy-plane after integrating over the z-axis. 

FIG. 6. Charge transfer upon the formation of the gold-BPD-gold contact. The top and the middle figure show the isosurface 
plot of the region where electron increases and decreases respectively. The bottom figure shows the spatial distribution of the 
transfered electrons as a function of position in the xy-plane after integrating over the z-axis. 

FIG. 7. Electrostatic potential change upon the formation of the gold-PDT-gold contact as a function of position in the 
xy-plane. Also shown is the projection of the molecule onto the xy-plane. 
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FIG. 8. Electrostatic potential change upon the formation of the gold-BPD-gold contact as a function of position in the 
xy-plane. Also shown is the projection of the molecule onto the xy-plane. 

FIG. 9. Band lineup at the gold-PDT-gold contact. The curves corresponding to the spin-up and spin-down electrons are 
virtually identical. The left figure plots the transmission versus energy, while the right figure plots the projected density of 
states onto the five frontier molecular states which can be identified from their peak positions. For PDT, these are -2.25(eV) 
(LUMO+1), -2.45(eV) (LUMO), -6.5(eV) (HOMO), -7.25(eV) (HOMO-1) and -7.7(eV) (HOMO-2). The horizontal line corre- 
sponds to the energy levels of four frontier molecular orbitals as plotted in Fig. 3. Note that the sharp peaks in the PDOS plot 
have been truncated here because showing them in full would have made the peaks corresponding to HOMO and LUMO less 
visible. 

FIG. 10. Band lineup at the gold-BPD-gold contact. The curves corresponding to the spin-up and spin-down electrons are 
virtually identical. The left figure plots the transmission versus energy, while the right figure plots the projected density of 
states onto the five frontier molecular states which can be identified from their peak positions. For BPD, these are -2.5(eV) 
(LUMO+1), -2.85(eV) (LUMO), -6.35(eV) (HOMO), -6.95(eV) (HOMO-1) and -7.45(eV) (HOMO-2). The horizontal line 
corresponds to the energy levels of four frontier molecular orbitals as plotted in Fig. 4. 



FIG. 11. Characteristics of the surface perturbed HOMO and LUMO molecular states at the gold-PDT-gold contact. 



FIG. 12. Characteristics of the surface perturbed HOMO and LUMO molecular states at the gold-BPD-gold contact. 

FIG. 13. Characteristics of the metal induced gap states. Left figure: gold-PDT-gold contact. Right figure: gold-BPD-gold 
contact. 

FIG. 14. I-V (upper figure)and G-V (lower figure) characteristics of the gold-PDT-gold device. The inset in the I-V plot 
gives the maginified view at low bias. The dotted line in the G-V plot is obtained assuming the transmission-energy relation 
to be bias- independent. 



FIG. 15. I-V (upper figure) and G-V (lower figure) characteristics of the gold-BPD-gold device as in Fig. 14. 

FIG. 16. Bias-induced modification of molecular levels at gold-PDT-gold junction (left figure) and gold-BPD-gold junction 
(right figure). We have also shown the position of the equlibrium Fermi- level Ef and the electrochemical potential of the two 
electrode Hl(r) in the plot. 

FIG. 17. Nonequilibrium occupation of molecular orbitals as a function of voltage for the PDT (left figure) and BPD (right 
figure) molecules in the molecular junction. Here we show the electron occupation of the HOMO-1, HOMO, LUMO and 
LUMO+1 states. 

FIG. 18. Spatial distribution of charge transfer and potential drop at the gold-PDT-gold and gold-BPD-gold contacts for 
bias voltage of 3.0(V). 



FIG. 19. Bias-induced modification of molecular states and transmission coefficent at voltages of 1.6(V), 3.0(V) and 3.8(V) 
for the gold-PDT-gold contact. The sharp peaks in the PDOS plot are not shown in full here. 

FIG. 20. Characteristics of field-induced modification of molecular states at V = 3.8(V) for the gold-PDT-gold junction. 
Left figure (right figure) shows the LDOS at energy corresponding to the peak position in the projected DOS of the LUMO 
(HOMO) states. 
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FIG. 21. Bias- induced modification of molecular states and transmission coefficent at voltages of 1.4(V), 3.0(F) and 4.0(V) 
for the gold-BPD-gold contact. The sharp peaks in the PDOS plot are not shown in full here. 
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